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The event mixing method is analyzed for the study of the event-by-event K/n ratio distribution. 
It is shown that there exists some correlation between the kaon and pion multiplicities in the 
mixed events. The K/n ratio distributions in the mixed events for different sets of real events 
^ ' are shown. The dependence of the distributions on the mean K/n ratio, mean and variance of 

I multiplicity distribution in the real events is investigated systematically. The effect of imperfect 

, particle identification on the K/tt ratio distribution in the mixed event is also considered. 
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I. INTRODUCTION 



The strangeness production in high energy heavy ion coUisions is an important topic. Because there is no strangeness 
content in the initial states of the colhsions, the strangeness in the final states is a result of the complex interactions 
J>. in the collisions. It has been proved that the strangeness in the final states is sensitive to the dynamics of the process. 
, The enhanced strangeness production in ultra-rclativistic heavy ion collisions may be caused by: (i) equilibrated gluon 
rich plasma phase [Q; (ii) baryon-junctions (iii) diquark breaking and sea-diquark (iv) strong color fields Q etc. 
' In particular, the enhanced production may be related to the formation of the long searching deconfined quark-gluon 
plasma (QGP) ||^ and can shed light on the time scale of chemical freeze-out and therefore carry information about 
the early stage of a heavy ion collision 1^ . 

In recent years, with the advent of large acceptance detectors, event-by-event (EBE) analysis [Q, ||, ^, |l^, |ll|, |lj, 
|l3t has been performed in many investigations. The philosophy of EBE physics is based on following speculation: 
Although the conditions for QGP production may be reached in every event, the fact that a phase transition from 
QGP state to hadronic final state is a critical phenomenon implies that deconfined QGP state may occur only in a 
Oh, very small sub-sample of events. So the fluctuations accompanying the phase transition will, in effect, be averaged out 
1^ ■ in the conventional ensemble analyses. The EBE analysis searches for fluctuations of observables at the event level, 
so it can be used to select interesting or anomalous event candidates jl^ with specific dynamical properties. So the 
EBE analysis can provide dynamical information which cannot be obtained from the traditional inclusive spectra. 

To estimate the amount of dynamical contributions to the distribution of EBE observables, the method of event 
mixing has been used to create reference events which may contain only statistical fluctuations and have the same 
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d ' inclusive (momentum and etc.) spectra as the original ones. To create mixed events one first sets up a track pool 
( "event- mixer" ) and inserts many (a few hundred or more) real events into the mixer. Then, by selecting particles 
from different original events in the pool one can get a non-real (mixed) event. In the process of constructing a 
mixed event only one track in any original event is allowed to be used. The number of particles in a mixed event is 
a random integer satisfying the multiplicity distribution determined from the original real events in the filled pool. 
Thus, the mixed events are random samples from the same inclusive track population as the real events, but there is 
no internal momentum correlations in them. So, the mixed events may be useful in the study of correlations due to 
some conservation laws, such as energy-momentum conservation law, since the constraint of such conservation laws 
plays no role in the mixed events. 

Recently the EBE analysis is used in experiments to analyze flavor (strangeness) fluctuations in high energy heavy 
ion coUisions [l^. In one of our earlier papers [ p^ we tried to calculate the trivial K/n ratio distribution 

for events having statistically independent multiplicity distributions for kaons and pions. For that pure statistical 
case, the K/t: ratio distribution was shown to be approximately Gaussian, and its width depends on the shape of the 
distributions of kaon and pion multiplicities. That result is the starting point for investigating any dynamical effect 
in the production of strange kaons, though it is not for the real situation in physical problems. In real experimental 
casCjthe multiplicity distributions for kaons and pions may be not independent, and the mixed events were used in 
[f6| , [r^ to estimate the statistical contributions to the ratio distribution. In calculating the EBE K/tt ratio, the 
only things relevant are the numbers of kaons and pions in a real and/or mixed events. Unlike in the study of mean 
transverse momentum fluctuations, where the inclusive transverse momentum distributions for real and mixed events 
are the same, the inclusive distributions of kaon and pion multiplicities are not necessarily the same for real and mixed 



events. One can see this point clearly in Sect. [II. To the best of our knowledge, it is not known yet what will be 



the K/tt ratio distribution in the mixed events before event mixing is actually performed when the ratio distribution 
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in real events is given. Such knowledge is important for us to understand the true effect of event mixing, ft will also 
speed up the real data analysis and save a lot of computing resources. We should investigate theoretically on which 
quantities in the real events the K/t: ratio distribution in the mixed events depends. The first question one should 
ask is: What is the dependence of the ratio distribution in the mixed events on that in the original events? If such 
dependence is strong, the mixed events can hardly be used to estimate the dynamical contribution in the real data. 

In this paper, we try to study systematically the effect of event mixing on the EBE K/tt ratio distribution. By K/t: 
we denote the ratio of positively charged kaons to pions in the events. For the sake of simplicity, we assume that there 
are only pions and kaons in the phase space under study. This paper is organized as following. After some general 
discussion in Sect. II, we will start in Sect. Ill from the simplest case by assuming that the multiplicity and K/t: 



ratio are both fixed in the original events. For this case, the effect of event mixing on the K/n ratio distribution can 
be treated with rigorously and analytically. Then we will study the influence of K/tt ratio distribution in real events 
on the event mixing in Sect. There it is assumed that all the real events have common fixed multiplicity but 
different EBE K/n ratios. SectT^ will devote to the event mixing for original events having a non-trivial multiplicity 
distribution but with a common fixed EBE K/t: ratio. Sect. VI is for the most general case for mixing of events with 
distributed multiplicity and distributed Kj-K ratio. In Sect. VII we will consider the influence of imperfect particle 
identification on the ratio distributions in the real and mixed events. Sect. VIII is for discussion and conclusion. 



II. GENERAL DISCUSSIONS 



We first give some general formula for the calculation of some quantities characterizing the Kj-K ratio distribution. 
Suppose the probability for an event (real or mixed) to have N — k pions and k kaons is P{N, k). Define a generating 
function 



k TT 



then we have a global K/t: ratio Rg as 



Rg ^ - G'^(l, I)/ G' {1,1) , 



(1) 



(2) 



where the subscript together with a prime represents derivative of the generating function with respect to the variables. 
In the expression, (• • •) represents average over all possible configuration. 

In an EBE analysis, one measures EBE ratio r^ = fc/vr for every event and observes the distribution of the EBE 
ratio. For such an analysis, one can get 
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then one can get the width of the EBE K/tt ratio distribution as ct = -^/(fc^/vr^) — (fc/Tr)^. In deriving last equation 
condition 1 — P{N, N) ~ 1 is used, which is quite a good approximation for current experimental analysis. 



III. MIXING OF EVENTS WITH FIXED MULTIPLICITY AND FIXED K/n RATIO 

We begin to study the effect of event mixing from the simplest case. We assume that every real original event has 
the same fixed numbers of kaons and pions. Then the A'/tt ratio in the original events is fixed, so does the total 
multiplicity N. We denote in this paper such a case by FF (fixed multiplicity and fixed AT/tt ratio). No chemical 
fluctuations in the real events are assumed. 

Then, we consider the mixing of the original events. By construction, the mixed events have the same total 
multiplicity N as the original ones. When a mixed event is constructed, N particles are selected from the mixer, and 
one picks up one and only one particle randomly from tracks of a randomly chosen original event. Since the numbers 
of kaons and pions are assumed fixed in the original events, the selected particle from an original event can be a pion 
or a kaon with definite probability. Assume the multiplicity ratio of kaons to pions is R, the probability for the chosen 
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FIG. 1: Distribution of EBE K/n ratios in the mixed events from a set of original events with fixed muhiplicity and fixed kaon 
contents. 



particle to be a kaon is R/ (1 + i?), and that for a chosen particle to be a pion is 1/(1 + R). Then, the probability for 
a mixed event to have k kaons and tt ^ N — k pions is prescribed by binomial distribution 



PFF(iV,fc) 



R 



1 + R 
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(5) 



It is obvious that the multiplicity distributions for kaons and pions in the mixed events are determined by the total 
(kaons plus pions) multiplicity N together and the K/tt ratio R in real events. One sees that the inclusive kaon (or 
pion) multiplicity distribution in the mixed events is very different from that in original events. 

The EBE K/n ratio distribution in the mixed events together with its dependence on the multiplicity N and the 
original ratio R can be obtained from Monte Carlo simulation. To get a mixed event one can generate N random 
numbers. If a random number is less than i?/(l + R), we say that a kaon is produced. Otherwise, a pion is thought 
to be produced. If the total kaon number in a mixed event is k, one can get a K/t: ratio for the mixed event 
rg = k/{N — k). In Fig. |l| the normalized r,, distributions P{re) are shown for i?=0.15 and 0.2 as examples. For each 
case 5 million mixed events are generated. For each i?, four values of N , 200,250, 300 and 350, are taken. Gaussian 
fits to all the distributions are shown in the same figure. For fixed i?, one can see that the width of the distribution 
decreases with the increase of N . For the same iV, the width increases with the increase of R. 

The mean values and widths of the distributions can be obtained from a Gaussian fit to the histograms in Fig. 
^, and the fitted results are shown with error bars in Fig. ^ as functions of inverse root of multiplicity, \/\fN. It 
seems that finite particle number effect plays the only role in the mixing process for both R — 0.2 and 0.15 since 
a/R oc 1/VlV is quite a good approximation for the points shown in the figure. The slope for R = 0.15 is slightly 
larger than that for R = 0.20. The mean and the width of the distributions can also be calculated from Eq. (0). To 
use our general equations Eqs. (||) and (^), one needs to have the generating function for distribution Eq. The 
generating function for such a distribution is 



GFF{x,y) 



Rx + y 
l + R 



N 



(6) 



The calculated results for six multiplicities N are also shown in Fig. ^ as solid circles linked by curves. One can 
find that the mean values from two different ways agree with each other quite well. But there is slight difference 
between the widths from the two ways. The slight difference indicates that the EBE K/n ratio distribution in the 
mixed events is not a perfect Gaussian. For smaller R, the difference between the results from two ways is smaller. 
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FIG. 2: Mean values Re and relative widths a/Rs of the EBE K/tt ratio distributions in the mixed events as functions of the 
inverse root of the event multiplicity. The points with error bars are fitted results of Fig. ^ to Gaussian distributions, and 
points linked by curves are from theoretical calculation. 



IV. MIXING OF EVENTS WITH FIXED MULTIPLICITY AND DISTRIBUTED K/tv RATIOS 

As the second simple case for investigating effect of event mixing on K/n ratio distribution, we assume that the real 
events have fixed common multiplicity N but different strange contents (FD). Thus the real events have a non-trivial 
K/tt ratio distribution. 

Assume that the N particles in a mixed event are picked up from real events labelled ii, 12, • ■ • , itv and that the 
corresponding K/n ratio in those events are r^j^, r^^, ■ ■ ■ , ri^^, respectively. Not losing any generality, we assume that 
the particles picked from events Ji, Z2, • • ■ , ifc are kaons and others are pions. Since we have no interest in the order of 
the original events from which the kaons are picked up, the probability for such a configuration is then 

f N\ ri ra 5_ (-7-) 

V k J 1 + n 1 + rz " ' 1 + rfc 1 + Tfc+i " ' 1 + TAT ■ 

Let the mean K/t: ratio in the original real events be R, and rewrite as R + di. Since the mixed events are 
constructed by randomly picking up particles from randomly chosen original events, the final probability for a mixed 
event to have k kaons and N — k pions can be obtained from above expression by averaging it over all possible 
configurations. Since {Sf^^^) ~ for m=0, 1, 2, • • •, one can get 

PF^{N,k)^PpF{N,k)- [Pff{N -l,k-l)~PFF{NM+0{4) , (8) 

(1 + Ry 

where CTq — {6f) is the width squared for the ratio distribution in the real events. As stated above, the value of CTq 
is determined by the involved statistical and dynamical contributions in the strangeness production. Normally, CTq is 
quite small, therefore the contribution from higher order terms can be safely neglected. From last equation one can 
see that the width of the EBE K/n distribution will be slightly larger than in the case with /ia;ec? multiplicity and fixed 
K/n ratio. From this equation the mean and width of K/n ratio distribution can be calculated. To compare with 
the calculation, we can also perform a Monte Carlo simulation for the case. We assume that the distribution of K/n 
ratios in the original events is a Gaussian with a relative width ao/R = 0.2. With the same choices of R and N as in 
Fig. 1^, eight distributions of the EBE K/n ratios in 5 million mixed events are shown in Fig. |^. One can compare 
this figure with Fig. [l] and find that the corresponding distributions are almost the same. This is not surprising since 
the width cto used in the simulations shown in Fig. ^ is only 0.03 for R = 0.15 or 0.04 for R = 0.2. Thus the ctq term 
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FIG. 3: Distributions of the EBE K/n ratios in the mixed events from original ones with fixed muhiplicity and distributed 
ratios. The K/tt ratio distributions in original events are assumed to have a relative width (Jo/R = 0.2 for illustration. 

in Eq. has little effect on the final results. The mean values and relative widths of the distributions are also very 
close to those shown in Fig. ^ and will not be discussed any more. The insensitivity of the relative width a /Re (or 
the ratio distribution) in the mixed events to the detail of the ratio distribution of original real events is crucial for 
detecting dynamical contributions in experimental data analysis from the event mixing method. Because the original 
distribution is a mixture result from dynamical and statistical contributions and the dynamical correlations in the 
original events are moved in the mixed events, the difference between the two distributions from real and mixed events 
can be used to quantify the dynamical contribution to the distribution in the real events. But it should be pointed out 
that the dynamical contribution can make the width of the ratio distribution larger or smaller, as have been found in 
0_So, we think it better to use 

— ('^mixcd - CTo)^ to measure the dynamical contribution instead of cj^jxcd — '-^o 



V. MIXING OF EVENTS WITH DISTRIBUTED MULTIPLICITY AND FIXED K/tv RATIO 

In real experiments the event multiplicity is normally not fixed but has a non-trivial distribution. The multiplicity 
distribution will certainly affect the K/tt ratio distribution in the mixed events, because in Sect. Ill we have seen that 
the ratio distributions for mixed events with different multiplicities are different. To understand the event mixing for 
real experimental events, it is an important step to investigate the mixing of events with some multiplicity distribution. 

Suppose that the multiplicity distribution in the original events is M{N) and that all the real events have the 
same fixed K/tt ratio R. We denote this set of original events by DF. Since the mixed events have the same total 
multiplicity distribution as the original ones, the probability for a mixed event to have k kaon particles and N — k 
pions is then 

PBF{N,k) = MiN)PFFiN,k) . (9) 

From this probability one can get the inclusive kaon (pion) number distribution by summing over pion (kaon) 
number. It is clear that the joint probability Popik + tt, k) cannot in general be equal to the product of the inclusive 
pion and kaon number distributions. One can see this point from a simple fact that 



(fcTr) ^^^kTTPDF{k + TT,k) 7^ ^^A:PDF(fc + 7r,fc) x ^ ^ 7rPcF(fc + tt. A:) = (fc)(7r) 

fc TT k TT k T! 



6 




FIG. 4: Mean values Re and relative widths a/R^ of EBE K/n ratio distributions as functions of parameter r in the Negative 
Binomial multiplicity distribution for mixed events. 



This point can be directly checked even for the case FF with Eq. (||). From that equation one can get 

(kn) = N{k) - , (fc)(7r) = N{k) - {kf . 

They cannot be equal to one another for mixed events. That means there is some (non- dynamical) correlation between 
multiplicities of kaons and pions in the mixed events though the mixed events are constructed by randomly choosing 
one particle from each randomly chosen real event. The correlation is caused by the demand that the mixed events 
have the same total multiplicity distribution as the real events. 

Despite such an intrinsic correlation between multiplicities of kaons and pions in the mixed events, we can go further 
to analyze the K/t: ratio distribution in the those events. 

Let 



^z^M(7V) = /(z) , (10) 



N 



and using Eqs. ^ and the generating function for Pof{N, k) can be written as 



EE-VPz.wfc+-,fc) = /(f^) . (11) 

It is clear that the ratio distribution in the mixed events depends on the shape of total multiplicity distribution 
f{z) and the mean ratio R. In many experiments the multiplicity distribution can be parameterized as a negative 
binomial (NB) or modified negative binomial (MNB) distributions, whose generating functions are, respectively, 

GMz)^{l-r{z-l))-^ 
l-A(^-l) 



Gmnb(2) 



l-r(z-l) 



The NB distribution can be regarded as a special case of the MNB distribution with parameter A = 0. To see 
the dependence of the ratio distribution in the mixed events on the shape of total multiplicity distribution or on the 
parameters in the generating functions, we fix i? = 0.15 and the mean multiplicity {N) = 270.13, the same as in 
[p5|, |l6|. For the NB distribution, {N) = rn, therefore after (N) is fixed only one parameter is left which can be 
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FIG. 5: Mean values Re and relative widths a /Re of EBE if/yr ratio distributions as functions of parameter k in the modified 
negative binomial multiplicity distribution. Four variances w for the multiplicity distribution are chosen for illustration. 



used to specify the variance of total multiplicity distribution. In our calculation we choose r to change from 0.05 to 
0.89. The calculated mean value and relative width of the ratio distribution in the mixed events are shown in Fig. IJ. 
The mean value seems independent of the value of parameter r, while the width decreases with r slightly. Since the 
relative variance of the multiplicity distribution has a linear dependence on the parameter r in the NB distribution, 
the dependence of the ratio a/ Re on the relative variance of the multiplicity distribution is similar to that shown in 

i 

For the MNB distribution, there are three parameters k, r, and A. In our calculation, we fix i? = 0.15 and the 
mean multiplicity (A^) = K,{r — A) =270.13, the same as in the case for NB distribution, and we choose the relative 
variance of the multiplicity distribution w = ((iV^) - {NY)/{N) = r + A + 1 to be 1.0, 2.0, 2.2, and 2.5. Then the 
left parameter k is changed from 100 to 300. The results for the mean values and the relative widths of the ratio 
distribution are shown in Fig. |^. The mean Re is the same constant for all the choice of parameters, while the relative 
width a I Re is independent of k, and increases a little with the variance w of the multiplicity distribution. The w 
dependence of the relative width ajRe is a direct consequence of the multiplicity N dependence shown in Fig. ||. 
An important observation is that the calculated width g jRe for the mixed events is much smaller than 23% given in 

m 0- 



VI. MIXING OF EVENTS WITH DISTRIBUTED MULTIPLICITY AND DISTRIBUTED K/tt RATIO 

The last and the most general case for our investigation of the mixing of real events is that the real events have 
both the total multiplicity and K/tt ratio distributions (DD). For given total multiplicity distribution M{N) the 
corresponding joint probability for a mixed event to have k kaons and N — k pions is 

PDD(iV,fc) = M{N)PYD{N,k) . (12) 

Because of the result in Sect. |^ that the ratio distributions in the mixed events for cases FF and FD are almost 
the same, one has reason to expect that the results from Eq.([l^) will be almost the same as those from Eq. (||). So, 
this case will not be discussed any more. 
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FIG. 6: Mean values Re and relative widths a/ R of the intrinsic K/n ratio distribution as functions of the inverse root of 
multiplicity for events without perfect particle identification. No event mixing is done here. 



VII. EFFECT OF IMPERFECT IDENTIFICATION ON EVENT MIXING 



In high energy heavy ion colhsions, the identity of a particle is obtained from its charge and dE/dx information with 
the help of the Bethe-Bloch function. The rate of energy loss dE/dx depends on the charge, mass and momentum of 
the particle in a complex way, and different particles may have very close dE/dx values in some momentum range. 
Therefore, normally one cannot make definite conclusion about a particle's identity from such information. Instead 
one can only say, for example, that a track is a pion or a kaon with a probability less than 1. This means that a track 
which is thought to be a kaon may in fact be a pion, and vise versa. If one tries to get K/ir ratio with the Maximum- 
Likelihood Method from the dE/dx information, as done in [p^ , such uncertainty will play a role in the K/t: ratio 
distribution. Since the efficiency of identification depends on the momenta of particles, it is extremely difficult to 
model the identification process. In the following we only show that the identification efficiency has important effect 
on the observed K/t: ratio distribution. To do that, we assume that the identification efficiency is the same for all 
kaons (or pions). 

For simplicity, we assume that only kaons and pions exist in the final states. Considering the identification uncer- 
tainty of a track from its dE/dx information, we assume that the probability for a thought-to-be kaon to be indeed 
a kaon is pi, and that for a thought- to-be pion to be a real pion is P2- We assume that the values of pi (^2) is the 
same for all thought-to-be kaons (pions). For an event with K thought-to-be kaons and N — K thought-to-be pions, 
the apparent K/t: ratio is i? = K/{N — K). But the real numbers of kaons and pions may change from event to event 
even if N and K are the same for all events, and the probability P{N, k) of having k real kaons and N — k real pions 
can be calculated. In such a case, the K/t: ratio will have an intrinsic non-trivial distribution, although it is assumed 
that every original event has the same numbers of thought-to-be kaons and pions. It is not difficult to prove that the 
generating function of the intrinsic distribution P{N, k) is 

Ctairinsicl^'y) = (Pi2^+(1-Pi)y)^((l-P2)a;+P2?/)^"'^ ■ (13) 



From this expression, we know that the real global K/tt ratio should be 

R _ (fe) _ P1R+{1-P2) , . 

"^'-{N^k) (l-pi)^ + P2' ^ ^ 

which may be very different from the apparent ratio R = K/ {N — K). In the limiting case pi = P2 — 1/2, i.e., when 
the kaons and pions can't be distinguished, the global ratio Rg will be 1, no matter how large R is. Of course, the 
global ratio is equal to the apparent one if pi = p2 = 1 , as demanded physically. 
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FIG. 7: Mean values Re and relative widths a/R as functions of inverse root of the multiplicity for mixed events without perfect 
identification. 



We are particularly interested in the width a of the intrinsic K/n ratio distribution relative to the apparent mean 
ratio R, a/R. In the following we discuss the case with pi — P2 = P- More general cases can be studied in the 
same way. For definiteness we take in the following the apparent ratio R = 0.15, close to the K/tt ratios observed in 
current high energy heavy ion experiments. The mean and relative width of the intrinsic ratio distribution are shown 
in Fig. ^ as functions of inverse root of the event multiplicity. This relative width calculated from Eqs. (p^), 

and (^) increases with the 1/\/N. Though linear relation between a/R and 1/Vn can still be seen, the relation 
(T DC 1/ ^/N is not longer a good approximation. For poorer identification (smaller p) the relative width is larger for 
the same N. 

Then we study the effect of event mixing on the K/n ratio for the case with imperfect identification. As in Sect. 



Ill 



wc first fix the total multiplicity N and the apparent K/tt ratio R. Using Eqs. (^) and (|l3|), one can get the 



generating function for the joint distribution of real kaon and pion multiplicities in the mixed events as 



Rp+{l-p) , p + Rjl^p) ■ 

■X H z y 



1 + R 



1 + R 



N 



(15) 



and the resulting mean values and relative width of the K/ir ratio distribution are shown in Fig. |^ as functions of 
inverse root of the total multiplicity N for several identification efficiency p. The mean values are again constants as 
in the case without the event mixing, and the relative width decreases with but is larger than in Fig. |^ for the 
same Vl^. As in Fig. ^ the relation cr cx 1/ Vn is neither a very good approximation though the curves are almost 
linear. An important observation is that a/R depends now strongly on the identification efficiency p and can be about 
23% for p a little less than perfect value 1.0. 

Now we can study the mixing of events with a multiplicity distribution and imperfect particle identification. The 
generating function for the joint kaon and pion multiplicity distribution in the mixed events can now be formulated 
using Eqs. ( |l5|) and (|lo|). When the multiplicity distribution is assumed to be a modified negative binomial distribu- 
tion, the generating function for the joint kaon and pion multiplicity distribution in the mixed events with imperfect 
particle identification efhciency p is 



'-^mixcdV'''i y J 



1 - Al 



ru 



with 



[Rp + {1 - p)]x + [p + Rjl - p)]y 
1 + R 



(16) 



With this generating function we can perform theoretical calculation. We fix the mean multiplicity (N) = K(r— A) = 
270.13 as before. From Fig. | we have seen that the ratio distribution has extremely weak dependence on parameter 
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FIG. 8: Relative widths a /R as functions of w, which is a measure of the variance of multiplicity distribution, for mixed events 
with a modified negative binomial multiplicity distribution and without perfect particle identification. 

K. So that we fixed it to be 300. Then we can change the relative variance, w ~ {{N"^) ~ (^)^)/ {^) = 1 + r + A, of 
the total multiplicity distribution as a free parameter from 1.9 to 2.9. Since Rf. in this case is neither measurable nor 
meaningful, we focus only on the relative width a / R in the following. The dependence of the relative width of the 
ratio distribution is shown in Fig. |^. The relative width of the K /t: ratio distribution increases slowly with w but has 
a very strong dependence on the identification efficiency p. For poorer identification the relative width is larger. The 
value of the relative width can be close to 23% given in |jl^ for quite high identification efficiency. For other total 
multiplicity distributions, the mean and relative width of the EBE K /t: ration distribution in the mixed events can 
also be calculated in a similar way. 

VIII. DISCUSSION AND CONCLUSION 

From above sections one can know how to derive the EBE K/tz ratio distribution for the mixed events once the total 
multiplicity distribution and global K/n ratio in the real events are known. Then we can understand why the event 
mixing method can be used to compare with results on K /t: ratio distribution from real events. Contrary to naive 
expectation, there may exist correlation between multiplicities of kaons and pions in the mixed events. Although the 
mean EBE K/tt ratio in the mixed events is stable for all cases with perfect particle identification, the width shows a 
quite strong dependence on the total multiplicity distribution and the global K/ir ratio. But the ratio distribution in 
the mixed events has an extremely weak dependence on the shape of EBE K /t: ratio distribution in the real events. 
This property enables us to estimate the dynamical contributions to the ratio distribution in the real events. For 
events without perfect particle identification, the intrinsic uncertainty of particle identity makes the relative width 
of the K /t: ratio distribution in the mixed events larger and dependent on the multiplicity distribution in a more 
complicated way. Currently observed relative width of the EBE K/ n ratio distribution seems to be a combinational 
result from statistical fluctuations (mainly due to the finite total multiplicity) and the imperfect particle identification. 
More theoretical and experimental studies are needed. 

Acknowledgments 

This work was supported in part by NNSF in China. Part of the work was done when one of the authors (C.B.Y) 
was working in Max-Planck-Institut fiir Physik, Munich, Germany, and he would like to thank the Alexander von 



11 



Humboldt Foundation of Germany for the financial support granted to him. He would also like to thank Dr. P. Seyboth 
for the hospitality extended him during his stay in Munich. He also would like to thank Drs. St. Mrowczynski and 
M. Gazdzicki for reading the draft and valuable comments. 



[1] E. Shuryak, Phys. Rev. Lett. 68, 3270 (1992). 

[2] S.E. Vance and M. Gyulassy, Phys. Rev. Lett. 83, 1735 (1999). 

[3] A. Capella and C.A. Salgado, Phys. Rev. C 60, 054906 (1999). 

[4] M. Bleicher, W. Greiner, H. Stocker, and N. Xu, Phys. Rev. C 62, 061901 (2000), and references therein. 
[5] P. Koch, B. Miillcr and J. Rafclski, Phys. Rep. 142, 167 (1986). 

[6] L.P. Cserni, Introduction to Relativistic Heavy Ion Collisions, J. Wiley, New York, 1994. 

[7] M. Gazdzicki and St. Mrowczynski, Z. Phys. C 54, 127 (1992). 

[8] M. Gazdzicki, A. Leonidov and G. Roland, Eur. Phys. J. C 6, 365 (1999). 

[9] M.L. Cherry et al. (KLM Collab.), Acta Phys. Pol. B 29, 2129 (1998). 
[10] M. Bleicher, M. Belkacem, C. Ernst et al., Phys. Lett. B435, 9 (1998). 
[11] St. Mrowczynski, Phys. Lett. B 439, 6 (1998). 
[12] St. Mrowczynski, Phys. Lett. B 459, 679 (1999). 
[13] F. Liu, A. Tai, M. Gazdzicki et al., Eur. Phys. J. C 8, 649 (1999). 

[14] R. Stock, Page 507 in Proceedings of NATO Advanced Study Workshop on Hot Hadronic Matter: Theory and Experiment, 

Divonne-les-Bains, France, 27 Jun - 1 July 1994. 
[15] G. Roland, talk given in Hirschegg '97 Workshop on QCD Phase Transitions, Hirschegg, Germany, 13-18 Jan. 1997, NA49 

Note Number 116. 

[16] F. Sikler, talk given in Quark Matter '99, Torino, Italy 10-15 May 1999, NA49 Note Number 200. 

[17] S.V. Afanasiev et al., NA49 Collaboration, Phys. Rev. Lett 86, 1965 (2001); C. Roland, Ph. D thesis, Frankfurt University, 
1999. 

[18] C.B. Yang and X. Cai, Int. J. Mod. Phys. A 16, 1227 (2001). 



